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ABSTRACT 

Context. It has long been suggested that circumstellar disks surrounding young stars may be the signposts of planets, and still more 
since the recent discoveries of embedded substellar companions. The planet-disk interaction may create, according to models, large 
structures, gaps, rings or spirals, in the disk. In that sense, the Herbig star HD 142527 is particularly compelling as, its massive disk 
displays intriguing asymmetries that suggest the existence of a dynamical peturber of unknown nature. 

Aims. Our goal was to obtain deep thermal images of the close circumstellar environment of HD 142527 to re-image the reported 
close-in structures (cavity, spiral arms) of the disk and to search for stellar and substellar companions that could be connected to their 
presence. 

Methods. We obtained high contrast images with the NaCo adaptive optics system at the Very Large Telescope in L'-band. We applied 
different analysis strategies to probe the presence of extended structures or point-like sources using both classical PSF-subtraction 
and angular differential imaging. 

Results. The circumstellar environment of HD 142527 is revealed at an unprecedented spatial resolution down to the sub arcsecond 
level for the first time at 3.8 \im. Our images reveal important radial and azimuthal asymmetries which invalidate an elliptical shape for 
the disk as previously proposed. It rather suggests a bright inhomogeneous spiral arm plus various fainter spiral arms. We also confirm 
an inner cavity down to 30 AU and two important dips at position angles of and 135 deg. The detection performance in angular 
differential imaging enables the exploration of the planetary mass regime for projected physical separations as close as 40 AU. The 
use of our detection map together with Monte Carlo simulations sets stringent constraints on the presence of planetary mass, brown 
dwarf or stellar companions as a function of the semi-major axis. They severely constrain the presence of massive giant planets with 
semi-major axis beyond 50 AU, i.e. probably within the large disk's cavity that radially extends up to 145 AU or even further outside. 

Key words, stars : T Tauri, Herbig Ae/Be - stars : individual : HD 142527 - protoplanetary disks - instrumentation : adaptive optics 



1. Introduction 

Since the discovery of the first exoplanet dMavor & Ouelozl 
fl99l around the solar-like star 5 1 Pegasi, more than 760 have 
been detected so far using different observing techniques, from 
radial velocity measurements, photometric transit, microlens- 
ing to direct imaging. These discoveries reveal a great diver- 
sity of planetary systems in terms of physical properties (atmo- 
sphere and interior) and environment (single or multiple planets 
around single stars, in a circumstellar disk or even circumbi- 
nary planets). A new paradigm about the formation, structure 
and composition of planets is emerging, wider than what we 
have learned so far from the Solar System. Moreover, most of 
the known detected giant planets orbit within 5 AU owing to the 
detection biases of radial velocity and transit techniques. The 
statistical properties of this population of close-in planets, the 
period-mass distribution, the planet-metallicity correlation to- 
gether with the densities of Hot Jupiter interiors (Udrv & Santos! 
l2007t iMavor et alJl201 lh favor a formation mechanism by core 



* Based on observations collected at the European Organization for 
Astronomical Research in the Southern Hemisphere, Chile, ESO : run 
087.C-0299A. 



accre tion followed by gas capture (CA hereafter, Pollac ket al.l 
IT996T) . 

At wider orbits, the distribution of giant planets is less 
well known. However, the discoveries of planetary mass com- 
panions to young, n earby stars in direct imaging (AB P ic b, 
IChauvin et all 12001 RXJ 1609 b. lLafreniere et al. 20081), and 
further more around Foma lhaut dKalas et alJl2008h and HR8799 
dMarois et al.l 120081 l2O10h challenge the CA model. This sce- 
nario suffers from exceedingly long timescales and low disk sur- 
face density at large separations. Alternative scenarii could oper- 
ate at wide orbits, either comb ined with CA such as planet-planet 
scattering dCrida et alj|2009l) and outw ard planetary migration 
(see for instance Mordasini e t alJl2012h or alternative formation 
proce sses must be in voked like gravitational instability (GI here- 
after|c ameronl 19781) or stellar binary mechanisms. These forma- 
tion processes should lead to different physical and orbital dis- 
tribut ions of giant planets dBolevll2009t iDodson-Robinson et al.l 
2009). Hence, new discoveries are needed to understand planet 
formation at all orbits. 

The technique of direct imaging technique also offers the 
possibility to detect giant planets in circumstellar disks to di- 
rectly study the planet-disk interaction and even dynamically 
constrain the planets's mass dChiang et alJ l2009: Lagrange et al.l 
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I2012ah . Indirect signs of planetary formation such as clumps, 
gaps, holes and spiral arms (see a review in Papaloizou et alJ 
120071) have been probably already o bserved within s everal 
proto- planetary disks (e.g AB Aur. IFukagawa et al.l [2004; 
TCha. iHuelamo et al.l l20Tu I Andrews et al.1 1201 ll SAO 206462. 
iMutoet al.ll2012l) . 

In this context, using high contrast imaging, we have ob- 
served the close environment of the star HD 142527, known to 
host an extended proto-pla netary disk. HD 14252 7 is a Herbi g 
star lWaelkens et al.1 d 1996b . classified as a F6IIIe dHouklll97 " 
The H ipparcos measurement has been revisited by Ivan Leeuwe 
d2007l) leading to a parallax of n = 4.29 + 0.98 mas, i.e. d = 
233^3 pc. However, HD 142527 has be en identified as a member 
of the star-forming region Sco OB 2-2 (|A.cke & van den Anckerl 
120041) or Upper Centaurus Lupus (de Zeeu w et al.lll999l) hence 
the distance would be d = 145 + 15 pc in both cases. We choose 
to adopt this latter distance since we consider the identification 
of the membership more reliable. The age and the mass of the 
star are also matter of debate since Fukagawa et al. (2006, F06 
hereafter) derived an age of 2.2*5 Myr and M - 1.9 + 0.3 M 
from its stellar UV /optical luminositi es and isochr one (from 
iPalla & Stahlerlll999l) whereas IVerhoeff et al.l d201 lb found an 
age of ^ 5 Myrs and M = 2.2 + 0.3 M based on a new estima- 
tion o f the stellar lumin osity an d pre main sequence isochrones 
(from lSiess eTaTll2000l) . Finally. iPecaut & Mamaiekld2010l) dis- 
cussed the age of F-type members of the Sco-Cen OB associa- 
tion and estimated it to be in the range [5; 20] Myr. We therefore 
considered an age between 2 and 20 Myr for further analysis. 
The Table Q] summarizes the stellar prope rties of HD 142527 . 
This system, resolved by F06 in NIR and lFuiiwara et al.1 ((2006) 
in IR, presents a circumstellar disk quite unusual with unprece- 
dented infrared exce ss (Fi R /F* = 0.98 and F NIR /F* = 0.32, 
Dom inik et al.ll2003b compared to other Herbig stars. Evidence 
of grain growth (deduced from thermal and scattered light obser- 
vations) togeth er with the high l evel of crystallinity (MIDI obser- 
vation s, Ivan Boekel et al.ll2004l and NIR spectrum, iHonda et al.l 
2009) reveal that the outer disk is characterized by highly- 
processed dust. Moreover, F06 notes the presence of a spiral arm, 
a prominent cavity up to ^ 150 AU, two asymmetric facing arcs 
within the disk and measured a stellar offset of ^ 20 AU with 
respect to the center of the disk. Models have been computed 
to fit the infrared excess , the geometry of th e inner disk and the 
mineralogy of the dust (Verhoeff e t al.ll201 lb . They found a disk 
cavity between 30 and 130 AU separating an inner disk embed- 
ded in a halo and a massive outer disk extended to 200 AU. They 
conclude that on-going giant planet formation is suggested from 
the highly evolved status of the disk surrounding HD 142527. 

To get a more accurate view of the close-in structures (cav- 
ity, spiral arms) of the disk at unprecedented spatial resolution 
and to search for stellar and substellar companions that could 
be connected to their formation, we carried out VLT/NaCo ob- 
servations at L ' band using angular differential imaging (ADI 
hereafter). We present the observing strategy and the data re- 
duction in Section |2] the fine structures and the geometry of the 
cavity and of the outer disk in Section [3] finally the detection 
performance and the constraints on the presence of giant planets 
in Section |4] 



2. Observations and Data reduction 

2.1. Observing strategy and observing conditions 

HD 142527 was observed on July 29 th , 2011 under good see- 
ing conditions in the course of the thermal imaging survey of 



Table 1. HD 142527 properties 



Parameters 



RA (J2000) 
DEC (J200) 

Proper motion RA (mas/yr) 

Proper motion DEC (mas/yr) 

Spectral Type 

K (mag) 

Distance (pc) 

Age (Myr) 

Mass (M G ) 



+ 15 : 56 : 41.88 
-42 : 19 : 23.27 
-11.2 ±0.9 
-24.46 + 0.8 

F6IIIe 
4.98 ± 0.02 
145 ± 15 " 

2-20 
2.2 ±0.3" 



Notes. (a) Based on new estimation of the stellar luminosity in Verhoeff 
etal. 2011. 



56 young and dusty stars (Rameau et al. 2012, in prep). We 
used the VLT/NaCo dLenzen et alJ 120031: iRousset etaiTl2003l) 
high contrast adaptive optics system and imager in angular dif- 
ferential imaging mode (ADI. lMarois et aLll2006l) . The L'-band 
(Aq = 3.8 /urn, AA = 0.62 pari) was used with the L27 cam- 
era (platescale =* 27.1 mas/pixel). Twilight flats were obtained 
for calibration. A first set of short unsaturated exposures with a 
neutral density filter (ND-Long) was acquired to serve as cali- 
bration for the point-spread function (PSF) and for the relative 
photometry and astrometry calibrations. A second set of satu- 
rated exposures, without the neutral density filter, was obtained 
to image the environment of HD 142527 with a high dynamic 
range. Both sequences were composed of various offset posi- 
tions on the detector for background subtraction and bad or hot 
pixels removal. A windowing mode (512 x 514 pixels) was cho- 
sen to allow short exposure times (DIT) and favor the temporal 
sampling using the NaCo cube mode. The parallactic angle vari- 
ations (An = 44.4 deg) are reported on Table |2l together with 
the exposure time, the number of frames per cube (NDIT), the 
number of cubes (Nexp) and the observing conditions. The star 
HD 181296 (AO, K = 5.01 mag), that was used as reference for 
classical PSF-subtraction, was observed with the same observ- 
ing strategy during the same night. 



2.2. Data processing 

All images were reduced using the pipeline developed for ADI 
observations at th e Institut de Planetologie et d'Astrophysique 



observations at the Institut de Flanetologie et d Astr ophysique 
de Grenoble (see lLagrange et alJ 2010t IBonnefov et al.l 1201 lb 
iDelorme et aHl2012t IChauvin et al.ll2012l) . All data were homo- 



geneously flat-fielded and sky-subtracted ; bad and hot pixels 
were also cleaned from the images. The data were then recen- 
tered and bad quality images were removed from the cubes. The 
result was a master cube containing all individual cleaned and 
recentered frames, along with their parallactic angles. 

To probe the inner part of the environment of HD 142527, 
we used different PSF-subtraction strategies either adapted for 
the detection of faint extended structures or to search for point- 
like features. These four methods are described below: 

1 . Spatial filtering combined with simple derotation (refered as 
nADI). An azimuthal average of the radial profile is esti- 
mated and subtracted from each image of the master cube, 
followed by a derotation and mean-stacking. A further step 
consists in an additional median filtering (fnADI hereafter). 

2. Classical ADI (cADI). The median of all individual reduced 
images is subtracted from the master cube, followed by a 
derotation and mean-stacking of individual images. 
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Table 2. Observing Log of HD 142527 in L'-band with and without the neutral density (ND), as well as for the reference star 
HD 181296. 



Star 


Date 


UT-start/end 


DIT 


NDIT 


Nexp 


Tr-start/end" 


(Airmass)* 


(Co)" 


<ro>* 








(s) 






o 




(arcsec) 


(ms) 


HD 142527 (ND) 


07/29/2011 


23:49/23:52 


0.2 


80 


10 


-15.4/-12.9 


1.05 


0.76 


1.6 


HD 142527 


07/29/2011 


23:59/01:03 


0.2 


100 


129 


-7.9/36.5 


1.05 


0.55 


2.1 


HD 181296 


07/29/2011 


03:23/04:27 


0.2 


100 


124 


-5.5/23.3 


1.15 


1.34 


0.9 



Notes. <a) n refers to the parallactic angle at the start and the end of observations. (A) The airmass, the seeing Co and the coherence time To are 
estimated in real time by the AO system and averaged here over the observing sequence. 



3. The LOCI algorithm (Lafreni ere et alj | 20 07). The residuals 
are minimized for each image using linear combinations of 
all data to derive the PSF contribution at a given location in 
the image. For LOCI, we considered optimization regions of 
Na = 300 PSF cores, the radial-to-azimuthal width ratio g = 
1, the radial width Ar = 1 and 2 x FWHM, and a separation 
criteria of 1 .0 and 1 .5 x FWHM. A slightly revisited version 
of LOCI was applied with an estimation of the reference that 
does not take into account the contribu tion of the region of 
interest (referred as to mLOCI, see [Lagrange e t al. 2012b; 
iMilli et al.ll20Tl 

4. Finally, classical PSF-subtraction (e.g. iMawet et all 12009). 
The observations were not designed to perform such reduc- 
tion afterwards but we found that HD 181296 was observed 
in our survey (Rameau et al. 2012, in prep.) 2 hours later with 
the same parallactic angle conditions. We therefore used it 
as a reference (see table [2]fbr the observing conditions). To 
optimally subtract the PSF to HD1 42527 's mastercube, we 
first selected science and reference frames with similar par- 
allactic angle variations (covering a total of 26 deg). From 
each individual science frame, the closest reference one in 
terms of parallactic angle (within 0.3 deg) was subtracted, 
after recentering and flux rescaling. This step enables the op- 
timization the subtraction of the telescope+instrument static- 
speckles. The PSF-subtracted science frames are then dero- 
tated and mean-stacked. Different methods were applied to 
estimate the flux scaling factor between HD 142527 and HD 
181296 : 1) by applying the same scaling to all frames, 2) 
by determining an individual scaling factor for each individ- 
ual frame considering the total flux between the science and 
reference frames, 3/ using two independent part of the im- 
age(below 0.7 and above 1 .5 of radii, to exclude the bright- 
est part of the disk) to derive an individual scaling factor for 
each individual frame. The best result was validated accord- 
ing to visual check. 

These four PSF-subtraction methods allowed us to reach op- 
timal detection performance at all separations, identify biases 
and non-physical structures. Although ADI (as with high-pass 
filtering) processes are well adapted for point-like detections, 
they are not for axisymmetric disks with low inclination (in the 
present case i 30°) since the disks are heavily self-subtracted 
(Milli et al. 2012 submitted). However they still enable the detec- 
tion of non-axisymmetric structures and sharp edges in the disk. 
On the other hand, nADI, low pass filtering and PSF-reference 
subtraction reveal large scale features, both in terms of separa- 
tion and size. Careful checks have been done with the Airy ring 
patterns while interpreting the observed structures. The results 
of these four methods are discussed below in the context of the 
disk morphology and of the presence of close companions in the 
environment of HD 142527. 



3. The circumstellar disk around HD 142527 

3.1. Disk morphology 

Our results on the disk morphology of HD 142527 are mainly 
based on the classical PSF-reference subtraction method, more 
adapted to the conservation of faint axisymmetric structures with 
low inclination. The complexity of this method is the determina- 
tion of the flux scaling factor to avoid the creation of artifacts. 
After various tests, we ended up with the fact that there is not 
an optimal and unique scaling factor for all radii. We therefore 
tried to subtract the PSF-reference with a radial scaling factor 
at each 1 pixel-wide annuli or at each few pixels-wide annuli. 
The results were not convincing because of the asymmetric and 
high contrasted structures leading to over-subtracted regions. 
Unresolved diffuse material seems to highly contribute to the 
PSF-flux observed at L ' in addition to the observed structures 
further away (see details below). 

Nonetheless, it is evident from all images to probe fine highly 
contrasted structures. The Figure[TJand Figure|2]show the images 
of HD 142527 reduced using the classical PSF-reference sub- 
traction and fnADI, cADI, mLOCI procedures within a reduced 
FoV of 4.2" x 4.2". However, the photometry of the features 
remains sensitive to the choice of the scaling factor. The images 
reveal new complex, non-axisym metrical structures in addition 
to tho se previously reported (F06; iFuiiwara et al. I I2006T: lOhashl 
2008). We discuss the main features below. 

- Spiral arms. One of the main results from F06 reveals two 
elliptical facing arcs towards the East- West directions and a 
displacement of the star with respect to the center of a fitted 
ellipse on the disk. Our data do not confirm this architecture 
but rather suggest spiral shapes. The fnADI image (Figure|2] 
left) remains the most adapted to see the following features. 
From PA = 170 + 4 deg to PA = 264 ± 2 deg, the arm goes 
closer to the central star, from 1.17 + 0.04" to 0.95 + 0.02", 
with a higher bending at 226 deg. A remarkable feature ap- 
pears at PA = 265 + 1 deg : the spiral arm shifts closer to 
the star of 0.11 + 0.02". This represents a strong rupture to 
the axisymmetric structure suggested by F06. When going 
North, the arm moves away (up to 0.91 + 0.03"). We esti- 
mate the arm width of about 0.43 " for PA = 227 deg after 
PSF-reference subtraction with sf = 2.2 (without data bin- 
ning). 

Considering the East part of the disk (see Figure [2] right), 
the separation increases from North West (0.87 + 0.03 ") to 
South East (1.1 1+0.05"). This spiral arm is not homogenous, 
showing peaks and dips. fnADI, cADI and LOCI highlight 
this inhomogenity by revealing blobs. We also estimate the 
arm width for PA = 47 deg of about 0.38 ". 
Slighlty further out of the main western spiral arm, we detect 
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Fig. 1. HD142527 L' image with a 4.2" x 4.2" field of view reduced with the PSF-reference technique and a scaling factor of 2.2. 
This factor seems to subtract better the close-in flux without an over subtraction of the most external parts. North is up, East is left ; 
flux scale is linear. 



N 



.■ fc - 



1 orcscec 



Fig. 2. HD 142527 L' images with a 4.2" x 4.2 " field of view. North is up, East is left. Left : fnADI using a high-pass filtering. 
Middle : cADI. Right : mLOCI using a mask to exclude the disk. Note that the linear flux scale is different for each panel to reveal 
all structures. 



at least two spiral arms which are born out at (1.1 ± 0.1 ", PA 
= 242 + 3 deg) and (1,01 + 0.04", PA = 261+3 deg) and 
move away while going to North. These spirals show up very 
well in cADI and mLOCI reduced images, up to the 9<x con- 
fidence level within the LOCI ones. However, the intensity 
is not uniform along these arms and their width are narrow 
compared to the ones in the PSF-reference images. These ef- 
fects due to ADI processing are detailed in Milli et al. (2012, 
submit.). 

Further out from the star, we confirm the presence of a spiral 
arm between (^ 1.71", 243 deg) and (^ 1.88", 302 deg) 
as observed in K and H-band in F06. As mentioned in previ- 
ous papers at shorter wavelengths, we confirm that this arm 
is fainter with increasing wavelength. 

The brightness asymmetry between the West and East di- 
rections It clearly comes out from all images that the West 
side of the disk is brighter than the East one. This observa- 



tion confirms the previous one by F06. At the same sepa- 
ration, the Western arm is brighter than the Easter one up 
to ^ 2 mag/arcsec 2 (see for instance the SBD profile in 
Figur^Hl, bottom). 

The prominent cavity. We confirm this cavity already ob- 
served by F06. It is re-detected with an outer edge at 0.87 ± 
0.06", i.e. 126 AU at 145 pc, along PA = 227 deg. F06 
presented the cavity to be elliptical and later simulations took 
into account this geometry. However, it is clear that the outer 
edge of the cavity differs from an ellipse. In particular, it is 
much more pronounced into the South West and extends fur- 
ther up to 145 AU when going South, as a consequence of the 
spiral shape described earlier. Considering the East part, the 
edge is located at 0.73 + 0.02", PA = 64 deg, i.e. 106 AU. 
Two azimuthal dips. We confirm the presence of two big dips, 
detected at North (PA e [351 + 2; 9 ± 5] deg and South-East 
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(PA e [148 + 6; 164 + 2] deg). The presence of point-like 
sources in these two dips is discussed below (Section l4~2l . 



or by an asymmetry between the forward scattering branch and 
the backward scattering one. 



3.2. Surface brightness distribution 



Distance to the star (pix) 
50 i t 60_ 

Profile towards 47 

Fit towards 47 

4.7±0.5 




1000 

Distance to the star (mas) 



2UU0 



Fig. 3. Observed radial surface brightness distribution on 
HD 142527 L ' PSF-reference subtracted image with a scaling 
factor of 2.2 towards the extended part of the disk (47-227 deg). 
Dashed lines corresponds to fitted straight line between 0.9 " 
and 1.3" for PA = 47 deg and between 1.2" and 1.8" for 
PA = 227 deg. 



We extracted the surface brightness distribution (SBD) of the 
HD 142527's disk using the classical PSF-subtracted image. We 
derived the radial profile to evaluate the extension, the slopes 
of the outer disk and the amplitude of features observed in the 
images towards several PA in 1 pixel-wide annulus. The result- 
ing SBD toward the largest extension of the disk (i.e. PA =47 
& 227 deg) are shown on Figure [3] The slope was measured 
between 0.87" and 1.3 " for PA = 47 deg and between 1.2" 
and 1.8 " for PA = 227 deg. We found that the SBD follows a 
simple power low with a slope of ^ ~4.6. The mass of gas in 
the disk being poorly constrained (Ohashi 2008), if the circum- 
stellar disk around HD 142527 is gas-free disk, such an outer 
distribution would be indicative of the dust dynamics d ominated 
by radiation pressure effects (e.g. The bault & Wull2008l) . 

To investigate the brightness asymmetry between the 
West and East side o f the disk, we used the GraTeR code 
dAugereau et al.l fl999l) to simulate, as a preliminary step, the 
disk as a ring following a radial volumic distribution of grains 
: {{r/r )- 2ai » + (r / ro) -2 "^) -0 - 5 . We chose r = 145 AU, a oul = 
-a ln = -3.7Q, the disk height scale £ = 1 AU at ro, a gaussian 
vertical profile with y — 2, a linear disk flaring with [3 — 1 and 
an inclination of i = 30 deg. The disk width thus obtained is 
0.48 ", which is larger than the one measured for PA= 227 deg, 
but 0.36" if we assume i = 45 deg. We simulated the disk with 
different scattering properties, from g = (isotropic) to g — 0.5. 
An anisotropy with g = 0.25 could explained the flux ratio of 2 
with i = 30 deg and also g = 0.15 with i = 45 deg. However, 
such configuration could not explain the geometry asymmetry 
by itself. It might be created either by an offset of the disk center 



3.3. Role of a close-in binary companion? 

Bille r et al.l (1201 2|) have recently detected a close-in low-mass 
star c ompanion candidat e, using the NaCo/Sparse Aperture 
Mode dTuthill et al.ll2010l) . This candidate has a predicted mass 
of 0.1 - 0.4 M for a physical projected separation of 88 + 5mas 
(^ 13 AU). This candidate should rely at 3.3 pixels with a 
contrast of 5.2 mag with the primary in our NaCo ADI/L'-band 
images. However, we did not detect any point source or insight 
of this companion due to the very small inner working angle 
and thus numerous residual speckles. If confirmed with a sec- 
ond epoch observation, the companion could be responsible for 
the non-axisymmetric features of the circumstellar disk around 
HD 142527. Indeed, the secular perturbations from a central bi- 
nary are known to create transient, outward propagating spiral 
structures in a circumbi nary disc (see for instance the numer- 
ical investigations of I Artvmowicz & Lubowl 1 1994b iBatd [2000; 
Ide Val-Borro et al] 1201 lb . The extent, evolution and longevity 
of these spiral structures depend on several parameters, such as 
the binary mass ratio and orbital eccentricity but also the spatial 
structure and density of the gas component of the disc. The de- 
tailed study of their evolution in the specific case of HD 142527 
requires advanced numerical modeling that clearly exceeds the 
scope of the present paper and is left to future investigations. 

4. Additional companions around HD 142527 

4. 1 . Detected point-like sources 



N 








■ li 


it 


cc-1 — 




1 orcsec 


cc-2 — ► ' 



1 The distribution slopes have been chosen to match the SBD slopes 
measured on the data. 



Fig. 4. fnADI with median filtering image of HD 142527 with a 
1 1 " x 1 1 " FoV. The arrows indicate candidate companions ; cc-1 
and cc-2 are confirmed as background sources whereas the third 
one is likely an ADI artefact created by the disk. Note that the 
negative parts (black) close to the disk result from the azimuthal 
averaging. 



We looked for point-like sources in a complete FoV of 
19" x 19", obtained after considering all fields observed during 
the dithering pattern. We confirm the presence of two faint point 
sources (the companion candidates cc-1 and cc-2, see Figure©, 
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Table 3. Relative Astrometry and Photometry of the point-sources detected around HD 142527. 



Candidate 


UT Date 


Instrument 


Separation 


PA 


AL' 




Status 








(") 


(deg) 


(mag) 






cc-1 


04/06/2004 


Subaru/CIAO 


4.23 ± 0.02 


140.8 + 0.3 










07/29/2011 


VLT/NaCo 


4.168 ±0.008 


138.9 ± 0.2 


13.3 ± 


0.3 


Background 


cc-2 


04/06/2004 


Subaru/CIAO 


5.56 ± 0.02 


218.5 + 0.3 










07/29/2011 


VLT/NaCo 


5.371 +0.009 


219.9 + 0.2 


11.9 + 


0.2 


Background 




Distance to the star (mas) 
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Fig. 5. 2D contrast in L' at 5cr for point-like sources from LOCI 
reduction around HD 142527. 



already detected by F06. Combining the Subaru/CIAO dataset 
from 2006 and our NaCo observations, we confirm the back- 
ground nature of these two sources. Their relative astrometry and 
photometry derived at each epoch are reported in Table [3] 

At closer separations within 2", various point-like structures 
are detected in the cADI and LOCI images (see Figure[TJ). Owing 
to the low inclination of the disk, the ADI processing may create 
artifacts, either extended ones (holes) or point-like (in the case of 
abrupt changes in the shape of extended structures). The bright 
signal located at 1.3" and PA = 227.5 deg) is for example cre- 
ated by the elbow's shape of the south-western part of the disk. 
Consequently, we do not detect any unambiguous companion 
candidates within 2". 



4.2. Detection limits 

To compute the ADI detection performance with LOCI and 
cADI, the flux loss of the ADI processing was estimated by in- 
jecting fake planets of given separation and contrast in a cube 
with empty frames (and in the data master cube for LOCI). The 
full pipeline was run. By comparing the injected fluxes to the 
recovered ones, we could derive the flux losses as a function 
of the separation. The second step was to estimate the noise in 
the residual maps using a sliding box of 9 x 9 pixels. The 5<x 
detection limits were then derived by correcting the noise from 
the flux loss and applying a normalization using the unsaturated 
PSF. We visually checked and confirmed the detection limits by 
inserting fake planets with several expected contrasts at several 
separations. The 2D-contrast map obtained with LOCI is shown 
on Figure|5] We are sensitive to point-like sources with a contrast 
of AL' = 10.5 mag at 700 mas (A proj = 101 AU) and 12.5 mag at 




50 100 150 200 250 
Distance to the star (AU) 



300 



Fig. 6. Detection limits expressed in Jupiter mass towards the 
South direction with cADI (similar performance as LOCI). The 
aqua zone displays the explored masses between 2 and 20 Myr 
since the age of the star is not well established. The position of 
the disk as been added towards the same PA. 



1500 mas (similar performance was obtained with cADI). They 
correspond to planetary masses of 6.5 Mj up and 2.5 Mj up at 101 
and 218 AU, respecti vely, using the pre dictions of the COND03 
evolutionary models (Baraffe et alj|2003l) for an age of 5 Myrs. 
We also explored the performance towards the northern promi- 
nent dip of the disk and along the South direction within the large 
southern cavity. The ID predicted mass vs projected separation 
(only one case is shown because similar performance is obtained 
in both directions) towards the South is reported in Figure [6] 
Considering the large uncertainty on the age of the star, whereas 
the Herbig status and gas presence are evidence for young age, 
the explored masses between 2 and 20 Myr are plotted on the 
Figure. We do not detect any planets more massive than 5 Mj up 
in the northern disk's dip and 4 Mj up in the sourthern cavity. 
Careful considerations have to be taken since the detection lim- 
its might be overestimated close to the edges of the disk and 
to the spiral arm s due to the disk impact in ADI processes (see 
Mill i et alJl2012h . 

4.3. Detection probabilities 

To set constraints on the physical characteristics of unseen gi- 
ant planets or brown dwarfs that could be responsible for several 
disk structures, we used t he Multi-purpose Ex oplanet Simulation 
System (t he MESS code Bonavita et al.l2012h . Such analysis ini- 
tiated by dKasper et al.l 120071) has been extensively done up to 
now, especially to constrain the giant planet population at wide 
orbits. We first assumed a population of planets coplanar with the 
disk by fixing the inclination to 30 deg. The code creates a planet 
population using a uniform (mass and semi-major axis) grid. We 
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Fig. 7. Detection probability map of a faint companion around 
HD 142527 as a function of mass and semi-major axis. The con- 
straint of i = 30 deg has been set to orbital properties of the 
generated planets. This was calculated using the LOCI 2D de- 
tection limits within a 4.2" x 4.2" FoV. 

explored a parameter space of semi-major axis ranging between 
1 and 500 AU and of masses between 1 and 20 Mj up , resulting in 
a final grid of 100 x 100 points. For each point, the code has gen- 
erated N gen = 10000 orbits and computed the projected position 
of the planet at the tim of the observation. Using the 2D-contrast 
maps converted in terms of masses (using the COND model's 
predictions, Baraffe et al. 2003) and projected physical sepa- 
rations, the code finally checked the detectability of each sim- 
ulated planet and determined their detection probability given 
their masses and semi-major axis (see Figure[7] Our NaCo obser- 
vations set constraints for planetary mass-objects down to semi- 
majo axis a > 40 AU and do not completely explore the large 
cavity between the outer spiral structures of HD 142527 and the 
inner disk located at less than 30 AU. The detection probabili- 
ties unambiguously exclude the presence of massive planets with 
M > 9 M Jup and M > 5 M Jup for semi-major axis of a > 130 AU 
and 200 AU, respectively. To test the impact of inclination in our 
simulations, we ran simulations for i e [10,50] deg, consider- 
ing that a third component will not be necessarly coplanar if the 
binary companion detected by Biller et al. (2012) is confirmed. 
The detection probabilities do not change significantly our pre- 
vious results. 

5. Summary and concluding remarks 

Our observations using VLT/NaCo adaptive optics thermal 
and angular differential imaging reveal the disk surrounding 
HD 142527 at an unprecedented spatial resolution down to the 
subarcsecond level. The disk shows a vast and asymmetric cav- 
ity between 30 AU and 145 AU as well as two important dips 
towards North and South East. We also resolve a bright and 
inhomogeneous spiral arm towards the West direction, several 
fainter spiral arms and structured features rather than two el- 
liptical facing arcs as proposed in previous studies. The obser- 
vations in angular differential imaging allow us to derive un- 
precedented detection performance of point-sources down to 
planetary-mass at 40 AU from the star. Monte-Carlo simulations 
combined with our detection map set stringent constraints on the 
presence of substellar companions orbiting HD 142527. We ex- 
clude the presence of brown dwarfs and massive giant planets 



in the outer part of the cavity and the presence of giant planets 
down to 5 Mj up further out. 

The observed structures indicate that the disk is in an active 
phase of global evolution. We have already discussed the possi- 
bility that the binary companion candidate might be highly re- 
sponsible for the creation of these features. However, other pro- 
cesses can be invoked such as one or more unseen planets under 
formation in the disk. 

Gravitational instability is known to create global spiral 
modes and in some cases fragmentation. However, given that 
the gas mass, as well as the surface density had been poorly con- 
strained so far, it is not possible to know if the conditions are fa- 
vorable for gravitational instability. Type II migration of a high 
mass plan et can also be a mechani sm able to produce such struc- 
tures (e.g. lArm ita ge & R ice 2005) in a gaseous disk. Thus, once 
again, constraints on the gas properties will allow to exclude or 
not the type II migration. 

While this article was under review Casa ssus et alJd2012l) re- 
ported Gemini/NICI Ks and L' observations of the outer part of 
the HD 142527 circumstellar environment. Given similar perfor- 
mance of NICI and NaCo, it gives us the opportunity to validate 
the observed features and excluded them to be instrumental or 
analysis artifacts. They confirmed the presence of several struc- 
tures : the spiral arms, the asymmetric inner cavity, the two North 
and South South East dips and the inhomogeneity of the East 
side. However, we do not confirm in our images the presence 
of what they called a 'knot' along the North West arm. In their 
images, we also clearly see the offset of the main arm that we re- 
ported towards the West, even if they did not mentioned this pe- 
culiar morphology. Furthermore, they investigated the presence 
of a 10 My protoplanet at 90 AU through hydrodynamic simula- 
tions to reproduce the morphology of the inner cavity and of the 
outer disk. Given our detection performance, such a companion 
would have likely been detected thus ruling out the input param- 
eters of this carving planet. They argued themselves for addi- 
tional ingredients within the simulations to reproduce all struc- 
tures for which one single planet cannot explain. Nevertheless, 
their approach is very interesting and could be improved by tak- 
ing into account our detection limits to set planet characteristics 
and the inner binary candidate companion within FARGO simu- 
lations. 

From the present data, it is rather difficult to identify the 
cause of the non axisymmetric cavity and spiral arms due to the 
lack of knowledge on the gas properties. However, we consider 
that they may be due to at least one bound-embedded object. 
Further observations and important numerical investigations will 
be necessary to further test the origin of the highly structured 
disk of HD 142527. 
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